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Abstract 

We present examples for the calculation of the distribution and fragmentation functions 
using the representation in terms of non-local matrix elements of quark field operators. 
As specific examples, we use a simple spectator model to estimate the leading twist quark 
distribution functions and the fragmentation functions for a quark into a nucleon or a 
pion. 

1 Introduction 

From the theoretical point of view, the distribution functions are quark correlation functions 
involving matrix elements of bilocal combinations of quark fields In fact, using the formal- 
ism of light-cone quantization, we can give nice physical interpretations for the leading twist 
functions fi{x), gi{x) and hi{x). The first one represents the probability of finding a quark 
with (light-cone) momentum fraction x. The second one measures the probability of finding 
right-handed quarks with momentum fraction x minus the probability of finding left-handed 
quarks with the same momentum fraction. It is called a chirality distribution or (if we work in 
the infinite momentum frame) helicity distribution. In this talk we shall not consider the third 
one. A similar interpretation can be given for the fragmentation functions Di{z) and Gi{z). 

The stucture functions that appear in the cross sections for hard scattering processes are 
weighted sums of these (and other) distribution functions. However, due to the non-perturbative 
nature of QCD at low energies, we are not able to calculate from first principles the distribution 
functions nor the fragmentation functions of quarks into hadrons. 

We are then lead to build models for the quark correlation functions. In this work we use a 
spectator model with some particular vertices that describe the interaction between the struck 
quark, the spectactor and the hadron. Our purpose is to calculate the leading twist distribution 
functions and fragmentation functions for the nucleon and the pion. 



We are going to neglect the contributions of gluons. They lead to and log corrections 
in the structure functions, the latter usually being absorbed in a scale dependence of the 
distribution and fragmentation functions. Thus, no QCD corrections will be computed. In this 
approximation the structure functions, the distributions and the fragmentation functions do 
not depend on Q"^ (for fixed x = —q^/P^ oi z = Pj^/k~). 



2 Distribution Functions 

It is well known that the cross-section for an inclusive deep inelastic lepton-hadron scattering 
can be written in terms of the hadronic tensor 

w^.. = ^ J^e^'HP,S\[UO,M0)]\P,S), (1) 

where M is the mass of the hadron, is its four-momentum, is its spin, is the momentum 
carried by the virtual photon and is the electomagnetic current. Since the hadron wave 
function is not known, we cannot calculate exactly the matrix element in (|1|). However, Wj^^, 
admits the standard parametrization in terms of the dimensionless structure functions Fi, F2, 
Gi and G2 (see, for instance, ref. 0). 

The leading contribution in Wfj^^ can be calculated from the well known handbag diagram. 
This diagram is composed of a hard scattering part (that can be calculated with perturbative 
QCD) and a soft hadronic part which must be modelled, since it involves unknown physics. 
This is the quark correlation function 

^.,{k,P,S) = / ^ e^'-« {P,S\ ^,(0) I P,S) 

^ J2{P,S I V^,(0) I X) 0(P° ) 5 ({k - Pf - Ml) {X I V^,(0) I P, S), (2) 
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where {X} is a complete set of intermediate states. A similar definition can be given for the 
antiquark correlation function 

Contained in this correlation function, is a set of distribution functions defined by 



$[ri(x) = ^ J dk-dkr Tr($r) 
where F is any 4x4 matrix. The leading twist distributions are 



(3) 

k+=xP+ 



fi{x) = <t>^^^\x) = ^ J dk-dkr Tr($ 7+) (4) 
= $[^^T^](x) = ^ / dk'dkr Tr{<l> 7+75) (5) 



9i[x 



and the structure functions are given by 



F2{x)/x = 2Fi(x) =^Qj [$[^^](a;) + $'=l^^](x)] (6) 



/ 



2AG'i(a;) = ^ [<l>[^''^^](x) + <l>^[^^^^l(x)] , (7) 
/ 



A being the helicity of the hadron. Moreover, it can be shown that G2 enters in the cross section 
muhiphed by a factor 1/Q (it is a twist 3 structure function) and then we shall neglect it. 
The most general expression for $ is 0| 



A. 



^k,P,S) = MA^ + A^f + A^l^+^-^a^^'P^K + iA^k-Si^ + MA^^-f^ 

+ ^k-S ri5 + ^k-S^j5+ tA.a^-'^.S.P, (8) 



+ z A,o a^" l5S^k, + z^k-S a^-" 75 k,P, + ^ e,,,^ YP^k^S''. 



Hermiticity requires that all the amplitudes Ai{k-P, k'^) are real. When time-reversal invariance 
applies, the amplitudes A^, Aq and Au vanish. 

The distributions can be written in terms of these amplitudes 

= TT J dadr e{xa - t - x^) {A2 + xA^) (9) 
gi{x) = n J dadr 9{xa — T — X ') [-^6- {Ar + xAs)^ (10) 

with T = k^/M'^, a = 2k- P/M^. 



3 Fragmentation Functions 

Also the calculation of the fragmentation functions involves a soft hadronic part given by 

A,,ik,P,,Sn)=Y: I tA e^'^ {Q\m)\Ph,SH,X) (X, P,, ^,|^,(0)|0). (11) 

This function admits an expansion similar to (^ but we are not going to need it. Instead, we 
define immediately the set of functions 



^^U) = — [ dk+d^PhT Tr(AT) 
Az J 



(12) 

k-=p-/z 



where z is the (light-cone) momentum fraction carried by the observed hadron originated by 
the fragmentation of the quark and PhT is the transverse momentum of this hadron in a frame 
where the hadron parent has no transverse momentum. 
The leading twist fragmentation functions are 



D,{z) ^ A[^-l(z) (13) 
Gi{z) = A^^'^^^\z) (14) 

and we can easily derive a useful relation between these fragmentation functions and the pre- 
viously defined distribution functions: 



D^{z) = zl2h{l/z) 
G^{z) = z/2g^{l/z). 



(15) 
(16) 



4 The Spectator Model 



In the spectator model the sum over all possible intermediate states in (H) is reduced to a single 
term, corresponding to the spectator \Xsp) . It seems reasonable to consider this spectator as 
a (scalar or axial vector) diquark. The amplitude for obtaining the state \Xsp), after removing 
a quark from the nucleon is given by 



{Xsp\ijm\p,s) 



u{p, s) 



It — m + ie 

We are going to assume that the vertices have the following structure: 



:i7) 



■pa 



V3 



Pa 

lal5- 75 



As form factors, we follow ref. 



and choose 



(18) 
(19) 

(20) 



Throughout this section, m stands for the (effective) quark mass, A is a parameter to adjust 
and A^ is a normalization constant. The other parameter of the model is the mass of the diquark. 
In order to fit the experimental data as well as possible, we set m = 0.36 GeV, m<j = 0.6 GeV, 
TUa = 0.8 GeV and A = 0.7 GeV. 

The choice a = 2 in (|20|) ensures the correct behaviour of the distributions at large x. For 
the case of a scalar spectator the correlation function is then 



U + m) V\k, P) if + M) (1 + 75^) 70 r^^(A;, P) 70 (^ + m)] (21) 

where is the mass of the scalar diquark. A similar expression is valid for the case of a vector 
diquark. 

Using equations (|1^) - (|21|) we can readily compute the amplitudes Ai of (^ and then 
calculate the distributions using (^ and (jT^) and the fragmentation functions using (jl^) and 



(|16|). A calculation of fragmentation functions for baryons has been given in Ref. 

For the pion one can also employ the spectator model. The spectator is now an antiquark. 
The correlation function can be written in terms of only four amplitudes (because the pion is 
a spinless particle). The vertex is now 

= ^ 75 (f + M) (22) 



with the same form factor as in (pOD but with a = 3/2 and A = 0.4 GeV. 

The preliminary results obtained with this model are shown in Fig. 3. At this moment no 
comparison with experiment is possible. This requires QCD evolution of these functions, what 
we postpone for a later work. At this stage the results ilustrate how the framework of quark 
correlation functions enables model calculations for distribution and fragmentation functions. 

This work is part of the research program of the foundation for Fundamental Research of 
Matter (FOM) and the National Organization for Scientific Research. 




Figure 1: Nucleon distributions and gi{x) (top); nucleon fragmentation functions Di{z) 

and Gi{z) (middle); pion distribution function fi{x) (bottom- left); pion fragmentation function 
Di{z) (bottom- right). For the nucleon, the full line is for the scalar diquark and the dashed 
line for the vector diquark. The normalization of the fragmentation functions is arbitrary. 
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